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We investigate the p-spin model with Gaussian-distributed random inter- 
actions in the microcanonical ensemble using the replica theory. For p = 2, 
there are only second-order phase transitions and we recover the results of 
Sherrington and Kirkpatrick obtained in the canonical ensemble. For p > 3, 
the transition between the ferromagnetic and paramagnetic phases is of first 
order, and the microcanonical and canonical ensembles give different results. 
We also discuss the ensemble inequivalence of the random energy model, cor- 
responding to the limit p — )• cxd. This is the first systematic treatment of spin 
glasses with long-range interactions in the microcanonical ensemble, which 
shows how the two ensembles give different results. 

1. Introduction 

Theoretical aspects of spin glasses attract interest not only from the statistical-physics 
community, but also from many other areas of science. The methods of spin glass 
theory have become indispensable tools for fields like statistical mechanics, information 
processing, image restoration and neural networks, see [HE] and references therein. Most 
models discussed in the mean-field context of spin glasses have long-range interactions, 
a typical case being the infinite-range Sherrington-Kirkpatrick model [3j. The purpose 
of the present paper is to investigate whether or not the canonical and microcanonical 
ensembles give different results on such systems with long-range interactions. 
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Interactions are said to be long range when they decay hke ~ 1/r", where r is the 
spatial distance between two interacting objects, and a is less than or equal to the 
dimension of the system, or when the range of the interactions is of the same order as 
the system size [3]. Systems with long-range interactions appear in virtually all areas 
of physics ranging from atomic physics to gravitational systems (see [H O [H O [8] for 
review) . Recent results on long-range interacting systems [H El [13 ttH [121 [13 151 [HI 
[3 [3 m [15] seem to indicate that we have to rethink the long-held principle that the 
physical properties of a system are independent of the ensemble they were derived from, 
especially where first order transitions are concerned [13 [13 [IB]- Related phenomena 
have been observed in driven systems with local dynamics |19t I2U]. The main reason for 
ensemble inequivalence in systems with long-range interactions is the lack of additivity 
of the energy: If two identical systems with energy E are brought together, the total 
energy is not IE because the interactions between the two subsystems are of the same 
order of magnitude as E. This suggests us to be cautious in the interpretation of the 
standard derivation of the canonical ensemble from the microcanonical ensemble. 

In a typical phase diagram of many materials, with the temperature and an external 
control parameter first-order phase boundaries are simple lines in the canonical 

ensemble, but in the microcanonical ensemble the question is sometimes more involved 
because there the temperature is not the principal control parameter. The temperature 
has to be derived from the equilibrium entropy, and this fact can sometimes have puzzling 
consequences, including negative specific heat, for systems with long-range interactions 
[2 H [22 l [23]. Experimentally, negative specific heat has been observed for a small cluster 
of atoms [Mj- 

Another feature of ensemble inequivalence is the appearance of a mixed phase instead 
of a clear phase separation in the region of a first-order transition [10 1 113 1 [T5]. Ergodicity 
breaking can also occur, when the thermodynamical phase space is not convex. This 
means that some of the intermediate states between two phase-space points are not 
realizable, inhibiting a smooth interpolation between these two points. Then, metastable 
states in a region of the phase space have diverging decay times since the system cannot 
smoothly evolve toward the true equilibrium state in another region of the phase space. 
When many-body interactions come into play, the ergodicity can already be broken 
before the phase transition occurs, inhibiting equilibrium phase transitions [23 [23115] • 

Spin systems on a lattice have been a tool of choice for detailed theoretical stud- 
ies since many models are tractable analytically and show many features of ensemble 
inequivalence [13 [13 [23 [23 [13 [2Z]- A paradigmatic, much studied system with long- 
range interactions is the infinite-range p-spin model, also known as the Ising model 
with infinite-range p-body interactions. It has been studied extensively in the canon- 
ical ensemble in the context of spin glasses. For, p = 2 the model is known as the 
Sherrington-Kirkpatrick model 0. For finite p, the p-spin model has been discussed 
in many papers, two of the most recent ones being [281 129] • while in the limit p — )■ oo 
it is known as the random energy model [3U] . In a recent paper [T5] we investigated 
the ferromagnetic p-spin model in random fields and have shown that the canonical and 
microcanonical ensembles lead to completely different phase diagrams for p > 3. To 
discuss ensemble inequivalence in systems with randomness in the interactions, includ- 
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ing spin glasses, however, new tools are required. In this paper we extend the replica 
trick to the microcanonical ensemble and apply it to the p-spin model with Gaussian- 
distributed random bonds and study whether the two ensembles lead to the same or 
different consequences. 

The paper is organized as follows. In section [2] we introduce the model and sketch 
the derivation of the replicated microcanonical entropy, with detailed calculations being 
given in the Appendix[Xl We compare the canonical and microcanonical phase diagrams 
for finite p and take the limit p — )• oo in section [31 We conclude with a discussion and 
summary in section HI 



2. Model and entropy 

In this section we introduce the model and explain the derivation of the microcanonical 
entropy and show its replica symmetric form. 

The p-spin model consists of N spins (in our case Ising spins Si = ±1, i = 1, 2, . . . , A^), 
where each spin interacts with all other spins through p-body interactions. Its Hamilto- 
nian is given by 

H = - ^ Jii,...,ipSi^ . . . Si^. (1) 

ii<...<ip 

In this paper the bonds (interactions) are quenched random numbers obeying a Gaussian 
distribution 

PiJn...^.) = (^) -P (- {Jn...^. - J^) ^) ' (2) 

with mean jo, called the ferromagnetic bias. For p = 2, this model is known as the 
Sherrington-Kirkpatrick model [3] and has been studied extensively in the canonical 
ensemble (see e.g. [T]). 

The thermodynamic potential in the microcanonical ensemble is the entropy, which 
is defined as the logarithm of the sum of states having some given energy, 5(e) = In $7, 
where we have chosen the unit fee = 1- Thus, we have to take the trace of a delta 
function over the spin configurations, 

n = Ti 6{Ne - H), 

where e is the energy per spin. Since we have to average over the quenched randomness 
of bonds, equation 1^, we employ the replica trick as detailed in Appendix |X1 For the 
entropy per spin, s = S/N, we get 

[fl"] = exp (nNs) , 

where we will have to take the limits A^ — )• oo and n — )> in the end. According to 
equation (fT3]) . the entropy per spin is given by 

ns = - ^(e + jo"iS)(Q"^)a/3(e + io"^p - ^ qai3qai3 " ^ maTfi^ + luTre^, (3) 

Q/3 a>l3 a 
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with 
and 

a>/3 j a i 

where the sums over a and /3 go from 1 to n. Here, = A^~^ 5'"S''' is the spin glass 
order parameter and nia = N~'^ - Sf is the ordinary magnetization. The symbols 
and fha arise from technical operations and can be understood as the Fourier modes of 
qafi and respectively. For details, see Appendix [XI 

We will work mainly with the replica-symmetric ansatz, = q and rua = m, and 
consider replica symmetry breaking where necessary. As shown in Appendix lA.H the 
replica symmetric entropy is 

SRS = -j2(^ + -^'o"^^)^ x -qP ~^ + j Dulnlcoshiy/^u + rh), 

where Du = du exp{—u'^ /2) / {Itt)^/"^ , with the saddle-point equations as given in equa- 
tions dSnii-d). 

To be able to compare the micro canonical phase diagram to the canonical one, we 
have to know the temperature of the microcanonical system. The inverse temperature 
is defined as the energy-derivative of the entropy, 

where m* maximizes and q* minimizes the entropy simultaneously. Notice that a pecu- 
liarity of the replica trick is that while the equilibrium entropy has to be maximal with 
respect to m, it has to be minimal with respect to q. 



3. Phase Diagrams 

In this section we compare the canonical and microcanonical phase diagrams of the 
model introduced in the previous section for finite p and in the limit p — t- 00. 



3.1. Sherrington-Kirkpatrick model 

The Sherrington-Kirkpatrick model {p = 2) [3j is well understood and has been studied 
extensively in the canonical ensemble. As only second-order phase transitions occur in 
this model, it is straightforward to confirm that the ensembles give equivalent results and 
lead to the same phase boundaries between the paramagnetic and ordered (ferromagnetic 
or spin glass) phases as long as the replica-symmetric solutions are concerned. 
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3.2. Many-body interactions 



The canonical phase diagram of the p-spin model with p = 3 is also well known [31^ 
I29j . The paramagnetic (PM) and ferromagnetic (FM) phases are separated by a first 
order transition. The spin-glass phase (SG) appears only when we consider replica- 
symmetry breaking (discussed in Appendix and the PM-SG transition is of second 
order. Between the FM and SG phases lies a mixed phase M where the FM solution 
shows replica symmetry breaking. The canonical phase diagram is depicted in figured] 
a). 



a) b) 




Figure 1: a) The canonical phase diagram of the model with p = 3. The equilibrium 
FM phase extends up to the blue dashed line and the spinodal line is drawn 
green dash-dotted, which marks the metastability limit of the FM phase. The 
horizontal thin dash-dotted line is the dynamical SG transition, b) The micro- 
canonical phase diagram. Here, the FM phase extends up to the blue dashed 
line, while the PM solution extends down to the red line, resulting in a mixed 
phase FM-I-PM. The green dash-dotted line is the microcanonical spinodal line, 
which limits the metastability of the microcanonical FM phase. Note that the 
microcanonical spinodal line (green dash-dotted) can lie below the equilibrium 
FM boundary (blue dashed), because the energy, not the temperature, is the 
control parameter in the microcanonical ensemble. The specific heat can be 
negative, thus a higher temperature may mean lower energy. For details, see 
text and figure [5J A dynamical SG-transition line as in a) may exist, but 
is omitted here due to lack of firm evidence. The black dashed line in both 
diagrams is the Nishimori line. 

Figure [1] b) shows the microcanonical phase diagram. The phase boundary between 
the PM and FM phases (red and blue dashed lines), as well as the spinodal line (green 
dash-dotted) were obtained by solving the replica symmetric saddle point equations 
()20b .-d). The PM-SG phase boundary (black) was calculated from the first-step replica 
symmetry breaking equations (|25h -e). In both ensembles the PM-SG transition occurs at 
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7pM-SG = 0.651. The SG phase extends up to the critical value of the ferromagnetic bias 
of j™^ = 0.768, where we have chosen the unit J = 1, and the metastability limit of the 
FM phase is at j™^ = 0.736. Between these two values of jo lies for low temperatures 
the region M with metastable FM solutions, while the SG and replica-symmetric FM 
phase meet directly around T = TpM-SG and jo = j™^ . Note that the replica-symmetry 
breakdown is used only for the SG solution in the present calculations. 

The microcanonical phase diagram deviates from the canonical one only along the 
PM-FM transition. The PM phase in the microcanonical ensemble extends down to 
the red line where the PM solution ceases to be stable (to be called the m = 0-line), 
while the FM phase extends up to the blue dashed line (m > 0-line), resulting in a 
FM-PM mixed phase (FM+PM). We can understand the reason for the emergence of 
the FM+PM phase in the microcanonical ensemble when we look at the entropy as 
a function of the energy as depicted in figure [2] a) . For each value of the bias jo there 
exists a FM solution to the saddle-point equations ([20b-d) which crosses the PM entropy 
at some critical energy. Since the inclines of both curves are different, this results in 
a temperature jump at the transition as shown in[2]b). Thus, the m = 0-line (red 
in figure [2 b)) corresponds to the temperature of the PM solution at the transition, 
while the m > 0-line (blue dashed in figure [T]b)) corresponds to the FM solution. The 
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Figure 2: a) Microcanonical entropy versus energy per spin for various values of jo) 

stands for the entropy of the paramagnetic solution, b) Corresponding caloric 
curves. The vertical dotted lines indicate the critical energies for jo = 0.8, 
jo = 0.9 and jo = 1.2. Negative specific heat occurs where dT/de < 0. 



appearance of a 'reentrant' tail of the FM transitions around the center of figure [T] b) 
implies an instability of the replica-symmetric solution. In the canonical ensemble this 
is well understood through the de Almeida-Thouless line [32] , but the generalization to 
the microcanonical ensemble is quite involved and will not be developed here, mainly 
because it has no direct relevance for the ensemble inequivalence. 

We can also draw the Nishimori line (black dashed in figure [T]) , which was recently 
generalized to the microcanonical ensemble |33] . The Nishimori line is a hypersurface 
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in the parameter space of models with quenched randomness, where the averages of 
physical quantities take on simple values [U |3l]. In the present case, the condition 
T = l/(2jo) defines the Nishimori line in the canonical ensemble and e = —jo in the 
microcanonical ensemble. An interesting feature is that on the Nishimori line there is no 
ensemble inequivalence [33] and the microcanonical m = 0- and m > 0-lines meet there. 

Part of the FM solution in the FM+PM phase for jo > 0.92 has negative specific heat 
closely below the critical energy. Below jo < 0.92 part of the metastable states have 
this property. To illustrate this point, we show the PM entropy and short segments of 
the FM entropy for various values of jo in figure [2] a) . The corresponding caloric curves 
are drawn in figure [2] b), where the metastable part is shown dashed, while the stable 
part is shown as a thicker line. Negative specific heat occurs when dT/de < 0. As jo 
decreases, so does the region of negative specific heat. Note that dT/de < occurs only 
in the mixed region FM+PM and never in the pure FM phase. 

The ergodicity is nowhere broken in this model. This is therefore an example for a spin 
system with infinite-range, many-body interactions in which not all features of ensemble 
inequivalence appear. 

The phase diagrams for p = 4 and p = 5 are qualitatively the same as for p = 3, while 
the region M becomes narrower and TpM-SG becomes lower as p grows. In tabled] we list 
the results of the critical values including results of the limit p — ?■ oo, discussed in the 
next section [3^31 





p = 2 


p = 3 


p = A 


p = 5 


p — )■ oo 


7pM-SG 


1 


0.651 


0.628 


0.615 


l/(2Vln2) = 0.601 


„-FMl 

Jo 


0.5 


0.726 


0.79 


0.811 


Vh^ = 0.833 


„-FM2 

Jo 


0.5 


0.767 


0.81 


0.82 


= 0.833 



Table 1: Summary of critical values for various p in the microcanonical ensemble. 



3.3. Random energy model 

In the limit p — )• oo, the model ([T]) is known as the random energy model [3D]. The energy 
levels are Gaussian-distributed random variables. It is also known to be the simplest 
spin glass [25] and its canonical phase diagram can be obtained analytically in a simple, 
straightforward manner. The SG phase in the canonical ensemble does not show up with 
the replica-symmetric ansatz, but is fully described by the first-step replica symmetry 
breaking (IRSB). The FM and PM phases are replica symmetric. 

Let us first derive the PM-FM phase boundary in the microcanonical ensemble. Here, 
we can still use the replica-symmetric ansatz. In the limit p — )• oo, the saddle-point 
equations (j20h -d) have the following solutions: q = m = (PM) and q = m = 1 
(FM), the latter being together with the requirement e = —jo. Taking the same limit, 
p — >• oo, for the IRSB saddle-point equations (p5b -e) it is straightforward to see, as in the 
canonical case, that the only other relevant solution, apart from the replica symmetric 
FM and PM solutions, is the SG solution gi = 1, ^o = and m = 0. To obtain the phase 
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diagram, we list the entropies of the different phases, 



spM(e) 

SFM(e) 



ln2-e^ 

(6 = 

1 



(- 

-io) 



^Ir2 < e < 0) 



ssG(e) = -(ln2-e") (-Vln2<e<0), 



(5) 
(6) 

(7) 



where x(0 < x < 1) is the IRSB parameter, the boundary between go and qi in the 
IRSB matrix (see Appendix IA.2I for details). In the energy range — Vln 2 < e < 0, the 
two entropies spm and ssg compete but the former wins because spm < ssG as long as 
< x < 1. As mentioned before, as far as the SG phase is concerned, the entropy has 
to be minimal. 




Figure 3: a) The canonical phase diagram of the random energy model with the energy 
e and the ferromagnetic bias jo as axes. The PM phase exists at energies above 
e > — \/ln 2 for jo < \/ln 2 and above the black line for jo > Vln 2. The SG 
phase exists at e = — \/ln2 for jo < Vln 2, shown green dash-dotted while the 
FM phase (red) has the energy e = —jo for jo > \/ln2. b) The micro canonical 
phase diagram. The PM phase exists for energies e > — \/ln 2, while the SG 
phase manifests itself only at the energy e = — Vln 2 (green dash-dotted) for all 
values of jo . The FM phase (red) shows the same behavior as in the canonical 
ensemble, it has e = —jo for jo > \/ln 2 and does not show up for a smaller 
bias, jo < Vln 2. The black dashed line in both diagrams is the Nishimori line, 
which coincides with the ferromagnetic solution (red) for jo > \/ln 2 



When X = 1, the two phases merge, but, on the other hand, x = 1 implies replica 
symmetry. Thus, there is no SG phase except at e = — \/ln2, where the two phases 
merge, since both have zero entropy. The FM phase is allowed only when the energy 
takes the specific value e = —jo . 

For any value of jo the FM and SG phases are defined only at a single energy re- 
spectively, and therefore it is not possible to define a temperature for these phases. We 
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cannot draw a standard phase diagram with the temperature and jo as axes to compare 
the microcanonical to the canonical phase diagram. We can, however, draw a canonical 
phase diagram with the energy and jo as the axes and thus compare the ensembles. This 
is done in figure [3l We see that the phase diagrams are almost equivalent, except that 
in the microcanonical case, figure [3] b), the PM extends down to e = — \/ln 2 for any jo. 

For a clearer understanding, we investigate how the microcanonical phase diagram 
evolves as p grows larger. To that end, we draw in figure H] the microcanonical phase 
boundaries between the PM and FM phases, scaled to their respective critical values 
(see tabled]) for p = 3 (black dash-dotted), p = 4 (red dashed) and p = 5 (blue). As the 
curves suggest, the area of the mixed phase FM+PM grows larger with p. We can infer 
that in the limit of very large p the pure FM phase extends only to the temperature 
where the PM phase freezes out for jo < j™^. The spinodal line and the m > 0-line 
coincide in the limit p — >• oo and the mixed phase M, of the FM and SG states, will 
disappear in this limit. The FM+PM phase extends over the whole quarter-plane right 
to the critical jo = \/ln 2 and above T = l/(2\/ln2). The SG phase is the same as in 
the canonical case, while the pure PM phase exists above T = l/(2-v/ln 2) and left of 
jo = \/ln2 as shown in figureOb). However, since neither the SG nor the FM phase have 
clearly defined temperatures when p = oo, the re-interpretation of figure [3] b) into the 
{T,jo) plane as drawn in figureOb) is to be taken with caution and has to be contrasted 
to results that were obtained in the microcanonical ensemble without the use of the 
replica trick 




Figure 4: The microcanonical curves show that the mixed region FM+PM, which lies 
between the two curves of the respective p, extends in an ever wider region for 
growing p. 



4. Conclusion 

We have derived the microcanonical entropy of the p-spin model with Gaussian-distributed 
random interactions with the help of the replica trick. To the best of our knowledge. 



9 



a) b) 




Figure 5: a) Canonical phase diagram of the random energy model, b) Microcanonical 
phase diagram as inferred from the limit of large p. The dashed line in both 
diagrams is the Nishimori line. 

our study represents the first consequent treatment of spin glasses in the microcanonical 
ensemble. For p = 2, we recover the results of Sherrington and Kirkpatrick, which was 
to be expected, as there are only second-order phase transitions present. For p > 2, 
many-body interactions effect a first-order transition between the PM and FM phases 
and part of the FM solution in the mixed phase FM+PM has negative specific heat 
above a certain value of jo, but never in the pure FM phase. This system does not show 
ergodicity breaking, which means that all metastable states have finite lifetimes for finite 
systems. We have also taken the limit p — )• oo within the framework of the replica the- 
ory and could show that the two ensembles give different results for the random energy 
model. 

Our study shows that great care has to be taken in the analysis of many-body spin 
glasses, and even in well-understood systems new properties can emerge. This fact 
should particularly be taken seriously in the spin-glass context, where the realistic finite- 
dimensional properties are often inferred from the analysis of mean-field limits. 
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A. Derivation of the Entropy 

In this section we will derive the entropy from the replicated number of states in detail 
and calculate the form of the entropy in the replica-symmetric and the first-step replica- 
symmetry breaking ansatz. 
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The Hamiltonian is given by 

H = — Ji^,,,i^Si-^...Sip, 



il<...<ip 

where the interactions are distributed via 



NP- 



-1 \ 1/2 



We want to calculate the microcanonical entropy per spin: s = \nO,/N where 

n = Tr5iE-H)=Tr [ ^^ME-h) ^ 

J 2ti 

and the trace is over the Ising spins. 

We use the replica trick and can closely follow the calculations for the canonical 
ensemble [Infi] = lim„^o ([^"] ~ 1)/^- The replicated sum of states reads 

m = 

I n dJ.,....Me-^^n...^.-^^)'^T^ (9) 

h<..<ip ^ \ a ^ J 

where a runs from 1 to n. We write (Yl^ dka/2TT) = Dk and integrate over the Jii.Ap to 
get 

[OT = Tr /z)fcexp[i^5^A;„-^^fc„fc^ ^ 5" ...^f^^,^ ...Sg ) 

\ a a, 13 h<...<ip J 



I \ 2 ]\^p-l ' 



(8) 



xexpli^j;^, E ^n-Sf, 



(10) 



a ii<...<i„ 



Now, we make use of the relation 



ii<...<ip \ i / 

which holds for ^ 1. We set J = 1 for simplicity and have 
[0"]=Trymexp|iVi^A:, ^i?/iV + jo ^ E )| 



fill 



X exp 



a,/3 
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It is convenient to introduce (5-functions — Sf) and (5(^^/3 — 1/^ S'f ) 

for a > /3. Then by rewriting = £^/A^ + jm^ and performing the integrals over k, we 
arrive at 



a i a>/3 i 

with Dq' = ]^ dg, £)m = H Q^^S = 1 for a = /9 and Qa/? = otherwise. Insert- 

ing the integral representations of the (5-functions and collecting the terms proportional 
to N allows us to evaluate the integral with the saddle point method as 

[J^"] oc expiV -^eoc{Q~^)afiep - X] 9a/9^a/9 -^rnarfia + InTre^ (13) 

y a/3 a>P a J 

with 

^ = E E sts^ + E E ■ (14) 

a>/3 i a i 

A.l. Replica-symmetric solution 

If we set = q and = m, we have 

Then the sum over the elements of the inverse of Qq,^ is calculated by a very simple 
trick: 



7 

E E('5~^)'y'9 = '^^ 



Noting that the sum of all columns of Qa^ yields the same result, J2/3 Qai3 = l+{n—l)qP, 
we arrive at 

ap 

We next have to trace over the Ising spins. Using 

L = i(T.^") +rn^S--^Q, (17) 

\ a / a 
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we can take the n — t- limit as 



If 1 

lim— lnTre^= / -DMln2cosh(-\/^M + ?fi) q, (18) 

n-^o n J 2 



with Du = exp(-«V2)/\/27f. 

The entropy defined by s = [r2"]/nA^ now becomes in the hmit n — )■ 0, — t- oo 

SRS = -j2(^ + -^o"^^)^ ]^ -qP ~^ \'^^~ + J Duln2cosh{y/qu + fh), (19) 
with the saddle point equations 

1~ J Du tanh^ {y/^u + fh) (20a) 

m = J Dutanh{^/^u + fh) (20b) 

Q = ii<^+hmn \!"^'~L (20c) 



,P' ' (1 -gP)2 



m = — 7o-(e + Jo^i^)-; r- (20d) 



A. 2. First-step Replica Symmetry Breaking 

The first-step replica symmetry breaking is characterized by splitting the diagonal of the 
matrix q^^^ into x x x blocks and filling the off-diagonal elements of these blocks with 
q\ and the rest of the matrix with q^^ while the diagonal is filled with 0. 

As such, the entropy is calculated from as s = ln[r2"]/nA^. The corresponding 
terms are calculated as 

n ^''^ ■'"^ 1 -I- (x - -I- (n - x)go 1 - (1 - x)(f^ - xq^ ' ^ ^ 



- (lapqap = nqoqo + x^{qiqi - qoqo)- - qiqi = -(1 - x)qiqi - xqoQo (22) 
n ^-^ X 

and 

- In Tr e^ = - / In / Dv cosh^ [uJ^^ vJqi - % + m) - -qi + In 2, (23) 
n X J J ^ '2 

where Dv is defined in the same way as Du. The entropy is 

(6 + jomP)2 , 1^^^ _ _ 

siRSB = ln2 - j2(-]^ - (1 - x)gf - x(f ) 2^^ a;)gigi - xgo^o - 9i) - 'mm 

+ - y Dnln y cosh^' (u\/^ -h v^Jqx -qQ^rn) (24) 
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and the saddle point equations are 



(a = 0,1) 



(25a) 



J2(i _ (1 _ ^)qP^ _ ^qV)2 



2pjo ^ ( e + Jo ) 



(25b) 



J2(l — (1 — x)q\ — xq^) 




(25d) 



(25e) 



(25c) 



K = u\/q^ + vy/qi - qo + m 
D'v = dv exp{-v^/2) cosh^ K. 



(25f) 
(25g) 
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